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Flow Equations for the BCS-BEC Crossover 
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The functional renormalisation group is used for the BCS-BEC crossover in gases of ultracold 
fermionic atoms. In a simple truncation, we see how universality and an effective theory with 
composite bosonic di-atom states emerge. We obtain a unified picture of the whole phase diagram. 
The flow reflects different effective physics at different scales. In the BEC limit as well as near the 
critical temperature, it describes an interacting bosonic theory. 
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Ultracold gases of fermionic atoms near a Feshbach 
resonance show a crossover [l[ between Bose-Einstein 
condensation (BEC) of molecules and BCS superfluidity. 
The controlled microphysics, which can be measured by 
two-body scattering and the molecular binding energy, 
and recent experimental breakthroughs can open a 
new field of quantitatively precise understanding of com- 
plex many body physics. On the theory side, this calls 
for a quantitative and reliable approach to strongly inter- 
acting systems. In turn, a precise experimental control 
of the relevant parameters, namely the scattering length 
a{B) depending on the magnetic field B, the density n 
and the temperature T, can test the viability of non- 
perturbative methods. 

The functional renormalisation group (FRG) directly 
connects the 'microphysics' to observable 'macrophysics' 
by a non-perturbative flow equation It has been 

used successfully for precision estimates in simple non- 
perturbative systems and has already been applied to 
coupled systems of fermions and collective bosonic de- 
grees of freedom in relativistic [J, Q and non-relativistic 
theories [1, 0|- In this approach, the results of perturba- 
tive renormalisation near the critical dimension Q or for 
a large number of components N Q can be recovered by 
an appropriate level of truncation of an exact functional 
differential equation. In a certain sense, the FRG can 
be regarded as a differential form of Schwinger-Dyson or 
gap equations in a 1PI L llj or 2PI setting, see (l2T |. 

Method and approximation scheme - We study the 
scale dependence of the average action [l3[. It in- 
cludes all quantum and thermal fluctuations with mo- 
menta q 2 > k 2 , or in the presence of a Fermi surface with 
effective chemical potential a > 0, all \q 2 — cr\ > k 2 . For 
k — ► 0, all fluctuations are included and T^o generates 
the 1PI correlation functions. In practice, this is realised 
by introducing suitable cutoff functions Rk(q) in the in- 
verse propagators. The dependence of on k obeys an 
exact flow equation 0], 



d k r k = isTr(rf 



R^dkRk 



(1) 



Here, STr sums over spatial momenta q and Matsubara 
frequencies lum as well as over internal indices and species 



of fields, with a minus sign for fermions. The second func- 

(2) 

tional derivative T k represents the full inverse propaga- 
te 

tor in the presence of the scale k. Both T k and T y k ' are 
functionals of the fields. 

In the present Letter, we demonstrate that already a 
very simple truncation of Pfe is sufficient to account for all 
qualitative features and limits of the crossover problem. 
We approximately solve Eq. (Q]) with the ansatz 



i^(d T - A - a)i' + ip*(d T - A v A)ip 
\-u((p) - h v (^p*il)-i.ij)2 - iftpttp^) ■ ( 2 ) 



r k = I d 4 x 



In addition to the fermionic fields i/j for the open-channel 
atoms, we use a collective bosonic di-atom field cp. De- 
pending on the region of the phase diagram and the 
scale k, it can be associated with microscopic molecules, 
Cooper pairs, effective macroscopic bound states or sim- 
ply represents an auxiliary field. The bosonic field is 
renormalised by a wave function renormalisation, (p = 

1/2 

Z v (p, such that at every scale k the term linear in the 
Euclidean time derivative d T has a standard normalisa- 
tion. (For the fermions, this renormalisation is omit- 
ted.) Eq. |T]) holds for fixed unrenormalised fields <p, i.e., 



( r L 2 iW = d 2 T k /dtp a d0p. We define 



o) 



(3) 



(2) 

where P^. is evaluated for an analytically continued 
Matsubara frequency wm — * ioj. The fields and cou- 
plings in Eq. ([2]) are scaled with powers of an appro- 
priate momentum scale k or energy scale k 2 /(2M) [Iol |. 
For nonzero density n, we choose the Fermi momentum 
k = fcp = {i'n: 2 n) 1 /' i . Our units are H = c = ks = 1- 

We consider a polynomial effective potential u(tp) writ- 
ten in terms of p = p*p>, 



m%p + \\p^ 
\\>{p- Po) 5 



SYM 
SSB 



(4) 



Here, we distinguish the symmetric regime (SYM), where 
the minimum of u is at p 



and rafp > 0, from the 



2 



regime with spontaneous breaking of the U(l) symme- 
try (SSB), where the potential minimum occurs at po(k). 
Superfluidity is signalled by po(k — > 0) > 0, with a gap 
for single fermionic atoms A = h^yfpQ. 

The flow starts at some microscopic scale k m with 



X^p — 0, m v in 



> and A v = 1/2. Here rr? ;„ is related 



to the magnetic field B and relative magnetic moment 
fi by dm 2 vin /dB = 2Mp/l b 2 , and reflects the detuning. 
We will concentrate on the limit of a broad Fcshbach res- 
onance, where /i 2 in — ► oo, m 2 in — > oo. In this limit, 
the microscopic action is strictly equivalent to a model 
containing only fermionic atoms with a point-like inter- 
action and scattering length a Then, the only rel- 
evant parameter is the concentration, c = afcp, (or ak 
for zero density), and the Feshbach resonance is located 
at a{B — ► Bq) — > oo. For broad resonances, the precise 
initial value of A v is unimportant. 

Finally, we specify the regulator functions for 
fermions and bosons. We work with optimised cutoffs 
fl2l. [l4| for space-like momenta (£ = q 2 — 0(a) a), 



Rl = Z v A v (2k 2 - q 2 )6(2k 2 ~ 
Rf = (fc 2 sgn£-^(fc 2 -|£D- 



(5) 



A central object is the flow of the effective potential u 
with t = lnfc/fcjn, displayed here for a < 0, 



fc 5 ( 7 

dtu = rj v pu' - —2 I — tanh7 v - 1 

07T \7 V 



(6) 



2^ 
3^ 2 v ' 



riAp+vA (ot+x 



■ coth a v — 1 



The functions 7, 7^, /3, a, a v , x read (for a < 0), 



7^ = v/y+T^, (7) 



fc 2 -q 
2T ' P ~ 2T 

2A v k 2 + u' pu" /—„ — 

— , ;>£=—, atp = y/ a* + 2xa. 



Primes denote derivatives with respect to p and the 
anomalous dimensions are rj v = —dthiZ v , tja^ = 
—dtliiAp. In our truncation, the Feshbach coupling 
/t, 2 = Z^h^ is independent of k. The flow equation is 
the analogue of similar equations in @. 

Vacuum limit - In order to make contact with exper- 
iment, we have to relate the microscopic parameters to 
the scattering length a for the two-atom scattering in 
vacuum. In our formalism, the vacuum correlation func- 
tions, that directly yield the cross section [l(|, are ob- 
tained from Ffe_>o m the limit n — ► 0, T — ► 0. For fixed 
fc the flow equations then simplify considerably. We find 
that for n = T = the crossover at finite density turns 
into a second-order phase transition [§, [l(| as a function 
of m 2 in or B, with 

m 2 > 0, <7a = atom phase (a~ J < 0) 

m 2 =0, oa < molecule phase (a^ 1 > 0) . (8) 

to?, = 0, (7a = resonance (a -1 = 0) 




FIG. 1: Chemical potential for T = minus half the bind- 
ing energy ?m/2 = — #(c -1 )c -2 . We compare our FRG result 
(solid) to extended mean field theory (dotted) and our previ- 
ous Schwinger-Dyson result (dot-dashed) [Io| . 



The dimensionless "vacuum chemical potential" cta = 
emM/k 2 is related to the binding energy eM of a molecule, 
see below. On the BCS side, the bosons experience a 
gap to 2 > and the low-density limit describes only 
fermionic atoms. On the BEC side, the situation is re- 
versed: fermion propagation is suppressed by a gap — <ta, 
and the low-density limit describes bound molecules. 
In the vacuum limit, we first solve the flow equation 

= 1), 



dtm 2 



K 



fc 5 



6tt 2 (fc 2 - a) 



(9) 



The condition that rh 2 vanishes for B = Bq, a = 0, k = 



leads to 



h 2 

2 _ - 2 _ V 7, 



2Mp 

k 2 



(B - B a ) - 2a. (10) 



In our picture, atom scattering in vacuum is mediated by 
the formation and decomposition of a collective boson. 
For the atom phase, one extracts the scattering length 
for k -> [H, 



h%k 



16nMfi(B - B a 



(11) 



Eq. (fTTj) relates /i 2 in = /i 2 to the scattering length a(B), 



thus fixing all parameters of our model. Eq. (fTTj) can also 
be used for B < B . Integrating Eq. ([9]) for a = a a < 
with the condition m 2 (k = 0) = yields the well-known 
relation between molecular binding energy and scattering 
length e M = cr A k 2 /AI = -l/(Ma 2 ). 

The flow of the renormalised Feshbach coupling /i 2 is 
determined by the anomalous dimension, 



h 2 



fc 5 



6tt 2 (fc 2 - a) 2 



(12) 



For a = 0, the rescaled renormalised Feshbach coupling 
rapidly approaches a fixed-point (scaling solution) given 
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A(T)/A(0) 



FIG. 2: Condensate fraction Qc (solid) and gap parameter 
A (dash-double-dotted) at T = 0. We compare Sic with ex- 
tended mean field theory (dotted) and Schwinger-Dyson equa- 
tions (dot-dashed) [10]. The condensate fraction matches a 
phenomenological Bogoliubov theory with bm = 0.92a in the 
BEC regime (dashed), consistent with our vacuum result. 



by %> = 1) h 2 /k = 6tt 2 . In the vacuum, we find A v = 
1/2 and Z v (a A < 0, k -> 0) = 1 + &*/(327iV = oa). 

Next, we study the equation for the dimensionlcss four- 
boson coupling \ v 



d t K 



k 5 



2V2XI A V {1 



T) V +T)A V 

5 



)fc 5 



4^ 2 (fc 2 -a) 4 3tt 2 (2Z v ,4 y fc 2 + m 2 ) 2 



(13) 

There are contributions from fermionic and bosonic vac- 
uum fluctuations, but no contribution from higher p 
derivatives of u. For a = and large /i 2 , we use the seal- 
ing form Z v = ^/(6vr 2 fc), m 2 = h%k = §, 

1]^ = 1, r)A v = and find for the ratio Q = X v k 3 //i 4 the 
flow equation 



The infrared stable fixed point Q* 
to a renormalised coupling 



. \ v _ 367r 4 Q, 
* ~ Z 2 ~ _ fc~~ 



(14) 

0.008 corresponds 
(15) 



<p 



to be compared with the effective four-fermion coupling 
A^, e ff = —h^/m 2 ^ = —6ir 2 /k. The constant ratio between 



these two quantities is the origin of the universal ratio 
between the scattering length for molecules and atoms, 
a>M/ a = 2A V / A^eg. 

In the molecule phase for a a < and k = 0, one 
has A^ jC ft = 8ir/y/— a a 0]. Omitting the molecule fluc- 
tuations, a direct integration of Eq. (fT3"j) yields A v = 
Sir/y/— a a and therefore am/ a = 2, whereas the molecule 
fluctuations lower this ratio. With the cut-off functions 
(|5"j) we get om/ci = 0.92, while further optimisation of 
leads to qm/o, = 0.71. Similar diagrammatic approaches 
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FIG. 3: Temperature dependence of the gap A(T)/A(0) in 
the BCS (solid, c _1 = -2), resonance (dotted, c 1 = 0) and 
BEC regime (dot-dashed, c _1 =4). 



give au/a = 0.75(4) [19(, whereas the solution of the 4- 
body Schrodinger equation yields om/b = 0.6 [l7| . con- 
firmed in QMC simulations flU ] and with diagrammatic 
techniques (la ]. 

Many-body problem - The system is now characterised 
by two additional scales, the temperature T and the 
Fermi momentum fcp. We set k = k-p from now on and 
use tildes instead of hats in order to indicate this spe- 
cific normalisation. We determine the initial values for 
the flow in these units by Eqs. (jTUJ) , (jTTJ) in terms of the 
concentration c = akp and h 2 , . For large h 2 v (broad Fes- 

hbach resonance), the value of /i 2 will not be relevant. 
Finally, we have to adjust er in order to obtain the cor- 
rect density, which is related to the a dependence of the 
potential at its minimum. Within our normalisation, this 
yields the condition du m { n /da = —l/(3n 2 ) for k = 0. We 
follow the flow of du m i n /da by taking the a derivative of 
Eq. ©, starting with an initial value — <7 3 / 2 6*(cr)/(37r 2 ) at 
k ln . The flow equation integrates out the modes around 
the Fermi surface for a > 0. At least for low T, the dif- 
ferent contributions on the right-hand side can be identi- 
fied with the densities in unbound atoms, molecules and 
the condensate [13] ■ Our result for <r(c _1 ) is shown in 
Fig. [TJ On resonance, we obtain cr(c _1 = 0) = 0.55, 
while quantum Monte Carlo (QMC) simulations give 
a{c- 1 = 0) = 0.44(1) [H, d{c- 1 = 0) = 0.42(2) 0. 

The density and temperature effects modify the flow 
when k ~ 1 or k ~ T 1 / 2 , i.e., when the wavelength of 
fluctuations being integrated out is comparable to the 
intcrparticle spacing or the de Broglie wavelength. For 
T = 0, in particular, to 2 reaches zero for fcsSB > 0, and 
the flow has to be continued in the SSB regime with 
p (k < /cssb) > until k — > 0. We show in Fig. [3] 
the condensate fraction Sic 03 and the gap for single 



po- In the BCS regime, the 



fermionic atoms A = / 
BCS value (A( C - 1 )/A sc ' s (c- 1 ) = 0.9) for the gap pa- 
rameter is approximately reproduced. On resonance, we 
find A(c~ 4 = 0) = 0.6, to be compared to the QMC value 
A^" 1 = 0) = 0.54 

At higher temperature, the effects of fermionic flue- 
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FIG. 4: Crossover phase diagram from the FRG approach 
(diamonds). We indicate the BCS (long dashed) and free BEC 
(horizontal line) values of T c and compare with Schwinger- 
Dyson equations (dashed and dot-dashed) [To|| . 

tuations on the build-up of po arc reduced and the 
bosonic fluctuations tend to diminish p . At T c where 
po(k — > 0) — > 0, we find a second-order phase transi- 
tion. The critical region is governed by boson fluctua- 
tions with universal properties in the 0(2) universality 
class. From the scaling solution, we find a critical expo- 
nent rj = r]tp + rjA v ~ 0.05 throughout the crossover. We 
plot A(T)/A(0) for different values of c" 1 in Fig.rj The 
universal behaviour is visible for T — ► T c . On the BCS 
side, the scale &ssb goes to zero for c^ 1 — > — oo, leading 
to an exponentially suppressed gap. 

The phase diagram in the (T, c _1 ) plane is shown in 
Fig. [D In the regime of weak attractive interactions, 
the BCS critical temperature is reproduced. On the 
BEC side, we find the shift of the critical temperature 
ATc/t bec = K ( nM )i/3 aM = (toP)-i/a K ( aM /a)c @ 
with k = 1.7, «m = 0.92a (short dotted line). Lattice 
simulations give k = 1.32(2) (2l| . 

In the BEC regime, both at zero temperature and close 
to T c the many-body physics reflects the behavior of "fun- 
damental" bosons of mass 2M interacting via a scatter- 
ing length om = 0.9a. This demonstrates the emergence 
of an effective bosonic theory, where all memory of the 
truly fundamental fermionic constituents has been lost, 
easily understood by the fact that the binding energy of 
the molecules is the largest scale in this region of the 
phase diagram. Moving to the unitary and BCS regimes, 



only a very narrow region around T c is dominated by 
bosonic fluctuations, which give rise to the above men- 
tioned critical behavior. Bosonic degrees of freedom are, 
however, crucial to accomodate the symmetry require- 
ments in the spontaneously broken phase throughout the 
crossover and form an important building block for our 
evaluation scheme. 

Our present truncation does not yet include the ef- 
fects of particle-hole fluctuations. They lead to a strong 
decrease of both the critical temperature and the gap pa- 
rameter in the regions of the phase diagram where there 
is a substantial Fermi surface, i.e., in the BCS and uni- 
tary regimes. In the BCS regime, the reduction of both 
T c and A(T = 0) bya factor (4e) 1/3 w 2.2 is well known 
as Gorkov's effect [22|]. In our formulation, this effect is 
encoded in the running of an effective four-fermion vertex 
which is generated by certain mixed boson-fcrmion dia- 
grams. The relevant diagrams precisely have the topol- 
ogy of the particle-hole diagrams in a purely fermionic 
setting. This will be discussed in future work. 

Conclusion - Our functional renormalization group 
analysis for ultracold fermionic atoms clearly demon- 
strates the necessity of the inclusion of bosonic quantum 
and statistical fluctuations beyond extended mean field 
theory. Both the BEC regime (c _1 — * oo) and the univer- 
sal critical behaviour (T — > T c ) are dominated by bosons. 
The vacuum fluctuations are crucial for the four-boson 
interaction. The thermal boson fluctuations are needed 
to establish the expected second-order phase transition. 
Our method is technically simple and involves only a few 
running couplings, still enough to resolve the full range 
of microscopic couplings, i.e., the BCS-BEC crossover, as 
well as the whole range of temperatures from the ground 
state to the phase transition. We control all regimes of 
densities including the physical vacuum (fcp — > 0) where 
the crossover terminates in a second-order vacuum phase 
transition. The simplicity of the picture constitutes an 
ideal starting point for systematic quantitative improve- 
ments by extending the truncation. For example, we have 
not yet included the (many-body) effect of particle-hole 
fluctuations which will lower T c in the BCS and crossover 
regimes. Extended truncations should lead to quantita- 
tive precision for the crossover physics. 

We thank S. Florchingcr, H.C. Krahl, M. Scherer and 
P. Strack for useful discussions. 
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